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Abstract 

We discuss the non-perturbative and the radiative corrections to inclusive 
B decays from the point of view known from QED corrections to high en- 
ergy e~^e~ processes. Here the leading contributions can be implemented 
through the so called "radiator function" which corresponds to the shape 
function known in heavy hadron decays. In this way some new insight into 
the origin of the shape function is obtained. As a byproduct, a parameter- 
ization of the radiatively corrected shape function is suggested which can 
be implemented in Monte Carlo studies of inclusive heavy hadron decays. 
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1 Introduction 

It has become a standard method to describe inclusive heavy hadron decays in 
the framework of an expansion in inverse powers of the heavy quark mass tuq. 
This approach allows us to compute lifetimes and branching ratios in a QCD 
framework, where radiative corrections due to gluon radiation can be systemat- 
ically computed, while non-perturbative effects can be parameterized to obtain 
estimates on their size. 

This method also yields reliable answers for differential decay rates such as 
energy spectra or invariant mass spectra that hold in practically all phase space 
with the exception of certain endpoint regions. These regions are characterized 
by the kinematic situation in which the heavy quark decays into a practically 
massless light quark, and all the emitted gluons are collinear to the light quark. 
In such a situation the hadronic invariant mass in the final state is small while 
the hadronic energy is still large. 

It is well known that in these regions the standard l/mq expansion breaks 
down and has to be replaced by an expansion in twists, very similar to what 
happens in deep inelastic scattering. The leading twist term has been discussed 
extensively |l[] and is described in terms of a single non-perturbative function, 
the light cone distribution function of the heavy quark. It has been shown that 
this approach is in for a certain range of parameters equivalent to the ACCMM 
model, which is very popular in the phenomenological description of inclusive 
heavy flavour decays. 

However, it is by not trivial to combine the radiative QCD corrections with 
this light cone distribution function in a consistent way. It has been observed in 
PI that one may think of the light cone distribution function as defined in |l[] as 
a convolution of a "soft" and a "jet-like" contribution, where the latter contains 
all the collinear singularities such as the Sudakov Logarithms. 

In the present note we give a practical parameterization of the radiatively 
corrected light cone distribution function which can be used to implement the 
leading radiative and non-perturbative corrections into a decay rate known at 
tree level. The spirit of this approach is very much along the lines which have 
been used in high energy e~^e~ or ep processes to implement the leading radiative 
corrections, e.g. for the Zq line shape 0. In these processes the the leading 
logarithmic resummation of the QED radiative corrections leads to a so called 
"radiator function" which is universal and has to be convoluted with the tree 
level process to obtain a radiatively corrected cross section. However, in the case 
of B decays the situation is a little more complicated due to the presence of non- 
perturbative effects and the Sudakov logarithms, but the physics is essentially 
the same. 

In the next section we discuss the formalism and introduce the parameteriza- 
tion for the analogue of the radiator function for heavy hadron decays. In section 
3 we consider a few applications of our formalism, finally we conclude. 



2 "Radiator Functions" in heavy hadron decays 

The situation in inclusive heavy hadron decays is completely analogous to the one 
in high energy e"'"e~-collisions. In both cases "initial state radiation" diminishes 
the cms energy of the "partonic" process. While in e~^e~' collisions the initial state 
radiation is the emission of soft and collinear (to the incident particles) photons, 
it is the emission of soft and collinear gluons in heavy hadron decays. The soft 
gluons correspond to a completely non-perturbative function, while the collinear 
ones lead to Sudakov-like Logarithms, for which also resummation formulae exist. 
In high energy e'^e^ processes the bulk of the QED radiative corrections is 
obtained by a convolution of the "partonic" (i.e. tree level) cross section with a 
radiator function, such that 

1 
d(j{s) = dxidx2 Dg-{xi,s,ml)Df,+ {x2,s,ml)datree{xiX2s) (1) 



where the tree-level cross section is taken at a reduced cms energy X1X2S. The 
radiator function D^-(+) (x) 

D{x, s,ml) = 6{l-x) + J ^^^P{x) + . . . ^ 6{l-x) + ^\n (^] P{x) + - ■ ■ 

ml 

(2) 
with 
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describes the energy loss of the incident particles due to radiation of soft and 
collinear photons including leading logarithms of the form ln(s/mg). The ellipses 
denote terms of higher order, which we shall not discuss here. Furthermore, the 
radiator function depends on the scale of the process ^Js and on a "regularization 
scale" nie cutting off the collinear singularities of the photon emission. 

While the full radiative corrections (i.e. the ones beyond the leading logs) 
have to be discussed at the level of amplitudes, the leading logarithmic radiative 
corrections are described on the level of probabilities: D{x) can be interpreted 
as the probability that the incident particles loose the energy fraction x due to 
radiation. 

It is well known that this concept is universal and can be applied to the QCD 
case as well. In deep inelastic scattering the radiator functions become the parton 
distributions and the scale corresponding to the lower limit of the integration in 
(0) becomes the factorization scale. 

In inclusive heavy hadron decays this works completely analogous. In this case 
the radiator function is the shape function or the light cone distribution function 
as given in 0. The physical picture is the same: The cms-energy available for 



(3) 



the decay is the mass M of the heavy hadron, but the partonic decay happens at 
a smaller energy niq + k^^ where now —niq < A;+ < A corresponds to the scaling 
variable x of the radiator function in (|^). 

The masses of the hadron M and of the mass uiq of the heavy quark are 
related to order l/mg by 

M = mg + A (4) 

The corresponding decay rates in the infinite mass limit are in analogy to ([l|) [|l| 

A 

dT= j dk+ R{k+, mg) dTtreeirriQ + k+) + 0(l/m^) (5) 

— oo 

where dTt^ee is the tree- level partonic rate and R{k+, fi) is the "radiator function" 
for heavy hadron decays. 

It has been argued in that this function can again be decomposed into a 
"soft" and a "jet-like" part, such that 

A 
R{k+,fi)= J dl+J{k+-l+,fi)f{l+,fi) (6) 
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where J can be computed in perturbation theory and / is a non-perturbative 
quantity. The factorization into a perturbative part and a non-perturbative con- 
tribution involves a factorization scale n, the dependence on which has been 
considered in [|, H. 

It is interesting to consider the dependence of this formalism on the on the 
heavy quark mass; we shall use this information later to adjust the parameters. 
It is well known that the total rates are proportional to the mass of the heavy 
quark ttiq such that 



r = 7^^ (l + b^^ + ■■■ + 0{l/ml)^ 



(7) 



where 6 is a constant depending on the process under consideration and also on 
the choice of the mass definition. In fact, using the pole mass in (|^ generically 
yields large and negative coefficients b and it has been argued that a different 
choice of the mass definition renders b small for all processes. For instance, using 
instead of the pole mass the MS mass and employing the one loop relation 
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the corresponding coefficient for the first-order correction to the rate becomes 
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On the other hand, we may also compute the total rate in the shape function 
formalism by integrating (^. Inserting the tree level result for the total rate, we 
have 

dk^ R{k+, ttiq) 7mM 1 + 5-^ + 0{l/ml)\ (10) 



We use the conditions on the shape function 
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dk+k+f{k+,ii) = Smqifi) 



(12) 



Note that the second relation is the equation of motion for the heavy quark, where 
a possible residual mass term drriQ 0, which represents the difference between 
the mass used in the Lagrangian and the mass used in the definition of the heavy 
quark momentum p = rnqv, has been kept. Obviously, at tree level we recover 
the total rate T = 'yrriQ, if we have Suiq = 

Our aim is to include the large radiative correction coming from the quark 
mass into a radiatively corrected shape function as defined in (^. This involves 
the first moment of R defined in (BI) and we obtain 
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(13) 
We assume that the non-perturbative function / satisfies the normalization con- 



ditions ([Tl|Jl^) at some small scale /i, and we find 
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/ rffc+J(A;+,/i) I (1 + 5^1+ f dk+^J{k+,fi) + 0{l/ml 



(14) 
This shows that the radiative corrections to the total rate are mainly given by 
the integral over the "jet-like" function. This is a universal function and hence 
we again find that some large portion of the radiative corrections (namely the 
leading log contribution) is universal. 

However, there also appears a contribution of order l/mq which according to 
(1^ has to be zero. This means that the residual mass term has to be chosen in 
such a way that the corrections of order l/mg vanish, i.e. 
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Let us now discuss the "jet-like" function explicitely up to order a^. The 
leading logarithmic piece is given by 
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J(A;+,/x) = (1 + ^6i)5(A;h 
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where the double logs \\i{—kj^/mq)/k^ are universal, while the single logs are not 
universal; their size is fixed by an appropriate choice of the scale yU. 

Using (|16|) we can compute the normalization and the first moment for which 

we get 
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where we have inserted —M as a lower cut-off for the ^+ integral. 

This cut off dependence reflects the fact that although all moments of the 
non-perturbative function / existQ, the moments of the radiatively corrected 
function R pick up power divergencies. The physical origin of this problem is 
the "bremsstrahlung" spectrum of the gluon emission behaving as l/{—k^) and 
thus the the flrst moment picks up a linear divergence. 

Clearly the physical limit of k+ is —mQ and thus we shall identify M = mq. 
The absence of l/mq corrections in the total rate thus means that the flrst 
moment of the non-perturbative function has to be 
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(19) 



while h\ is given by the large and universal correction to the total rate. 

Again we consider using a different mass deflnition in the shape function 
formalism. For instance, using the MS mass deflned in (^) we get for the mass 
difference 
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which at least for the non-logarithmic terms agrees with drriQ computed from 
the flrst moment of the shape function. This leads to the conjecture, that the 
bulk of the radiative corrections at the leading twist level can indeed be absorbed 
by an appropriate choice of the mass deflnition, even for differential quantities 
affected by the shape function. The shape function itself depends on the mass 



^This means that / has to decrease exponentiahy as k-^ 
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definition via its first moment, which has to be shifted by a Stuq corresponding 
to the difference between the mass definition employed and the pole mass. 

For practical applications, (p!6D is still not acceptable, since the "jet-like" 
function is interpreted in a similar way as in the case of e"'"e~ annihilation as a 
probability of having the tree level process at a reduced cms-energy. However, 
(p!6D is not positive definite. This problem has been discussed in ^ where it was 
shown that sub-leading contributions render the decay rates positive. Thus for 
phenomenological applications we add terms corresponding to sub-leading pieces. 
Of course, in general terms of sub-leading twist depend on the specific process, 
but here we add universal terms which are fixed by comparing to the results of 
p. The formula we are going to use is 



a. 



J{k^, ^) = 1 + ^61 5(fc^) + ^e(-A;+) 
on / on 
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The three constants appearing in (pl[ ) have to be determined from some input. 
One condition is that the large radiative corrections are taken into account, which 
fixes hi as discussed above. Furthermore, the minimal, kinematically allowed 
value of k+ is —mq where the rate should vanish, thus we require 



J{k+ = -mQ,fi) = 0. 



(22) 



Furthermore, the lower limit of the fc+ integration is now mq. 

These conditions can be fulfilled with suitable 61 and 63 for any choice of 62- 
The coefficient 62 multiplies a term of sub-leading twist that could in principle 
be neglected; for a vanishing 62? however, J{k^,fi) becomes negative for some 
—rriQ < k+ < 0. Comparison with explicit calculations ^ suggests a value for 62 
which in turn fixes the constant 63. This way we obtain 



J(A:+,/i) = (l-^b + 8lnJ^-4]6{k^) + ^ 
V Svr niQ J 3n 
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'k+ -3 ln(/i/mQ) + 2 
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Here b is from the 0{as) correction to the total decay rate (|^). 

The scale fi is not fixed by these considerations; explicit next-to-leading order 
calculations show that /i is about 0.1mQ...0.2mQ and thus is a low scale. Again 
this is similar to the high-energy e~^e~ case. 

To parameterize the non-perturbative contributions we use the model shape 
function from 1911, 



fik+) 



32 




n2 



exp 









(24) 



This function has vanishing first moment as required at tree level, i.e. Sniq = 0. 
A non-vanishing first moment as required once radiative corrections are taken 
into account is achieved by shifting the argument of / by fc+ — » fc+ + Smq, where 
delta rriQ is calculated from the first moment of (pT]) . 

The radiatively corrected shape function R{k+, /i) for a typical fi = 0.1146mfe 
(as obtained from an explicit calculation of the 0{as) radiative corrections to 
b —>■ uiui) is shown in figure |l|. Since for the decay B —>■ X^'-f the photon energy 
spectrum is (except for factors) just the shape function itself with the endpoint 
shifted to M/2, this diagram at the same time shows the spectrum for B -^ Xs^. 
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Figure 1: Radiatively corrected shape function R{k+,^). 



In the case of the decay B -^ Xu£ij the spectrum is a step function at maximal 
lepton energy, which is smoothened by the shape function effect. The correspond- 
ing results are shown in figure 



3 Applications to the decay B -^ Xgi^i 

As an application we employ the radiatively corrected shape function for the 
lepton energy spectra in i? — »• Xs£~^i~. While it is well known that the shape 
functions do not play a significant role in the leptonic invariant mass spectrum, 
they significantly modify the energy spectra of the leptons. 




Figure 2: Lepton energy spectrum for B -^ X^ii^, tree level (dotted line), and 
including corrections parameterized with the radiator function formalism (solid 
line). The rate is given in units of Gpiril \V^fj\ /9Qii^ 



Of particular interest are the lepton energy spectra, where the flavour of the 
decaying B meson has been tagged. These spectra have been suggested to test 
the coefficients Cg and Cio of the effective weak Hamiltonian governing weak 
decays of b quarks. In addition, for experimental reasons an energy cut has been 
applied to the second lepton. 

The partonic spectra show pronounced, step function like structures which 
will be affected strongly by the shape function effects. In figure (H) we show the 
Ei+ spectrum with a cut of 0.2mf, on the energy of the negatively charged lepton 
both with and without shape function effects. It can be seen how the radiatively 
corrected shape function not only shifts the endpoint of the spectrum and adjusts 
the rate but also smears out the features of the partonic spectrum. 

4 Conclusions 

The question of how to implement non-perturbative contributions due to bound 
state effects in combination with a consistent treatment of radiative corrections 
has been discussed form a new point of view. Comparing the situation in B decays 
to the case of e~^e~ annihilation we look at the shape function as the analogue 
of the radiator function in high energy e~^e~ processes. In both cases the bulk 
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Figure 3: Spectra for the energy of the positively charged lepton in i? ^ Xgd^d' 
with (sohd) and without (dashed) shape function corrections. Both spectra in- 
clude a cut of 0.2mb on the energy of the negatively charged lepton. The spectra 
are given in units of G']^a'^ \^ts^tb\ m^/3(27r)^ Standard values have been used for 
the Wilson coefficients; all our plots use ^^(mfe) = 0.215, and the shape function 
parameter is a /mi, = 0.57/4.71. 



of the radiative corrections (and in B decays the main non-perturbative effects) 
can be accounted for by a convolution of the tree level rate with an appropriate 
function, which can be interpreted as a probability for a certain energy loss. Since 
this works at the level of probabilities, this approach can be easily implemented 
into Monte Carlo programs, which has been done quite successfully for the case 
of high energy e+e^ processes 
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For the case of B decays the situation is complicated by the presence of non- 
perturbative effects. These have to be parameterized in terms of a so called shape 
function, which then is combined with the leading radiative corrections into the 
analogue of the radiator function. 

The energy available in the B decay process depends on the mass of the heavy 
quark in the initial state. On the other hand, the b quark mass is not a physical 
quantity and thus one has to consider the dependence on the mass definition. It 
turns out that the main part of the radiative corrections can indeed be absorbed 
into an appropriate mass definition; this fact is well known for the total rates, 
but it seems to be true also for the leading twist contributions parameterized by 



the shape function. 

We have given a parameterization of the radiatively corrected shape function 
which includes some sub-leading twist contributions which have been fitted in 
such a way that the explicitely known decays of the heavy-to-light type are well 
reproduced. This parameterization may be useful for practical applications, in 
particular for Monte Carlo simulations of B decays. 
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